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Abstract 

From the very early days of Particle Physics, both experimental and theoretical studies 
on proton-proton collisions had occupied the center-stage of attention for very simple and 
obvious reasons. And this intense interest seems now to be at peak value with the onset of 
the Large Hadron Collider (LHC)-studies at TeV ranges of energies. In this work, we have 
chosen to analyse the inclusive cross-sections, the rapidity density, the K/tt and p/7r-ratio 
behaviours and the < pr >- values, in the light of the Sequential Chain Model (SCM). And 
the limited successes of the model encourage us to take up further studies on several other 
aspects of topmost importance in particle physics with the same approach. 
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1 Introduction 

Proton-proton collisions are known to be the most elementary interactions and form the very 
basis of our knowledge about the nature of high energy collisions in general. Physicists, by 
and large, hold the view quite firmly that the perturbative quantum-chromodynamics (pQCD) 
provides a general framework for the studies on high energy particle-particle collisions pQ. Ob- 
viously, the unprecedented high energies attained at Large Hadron Collider (LHC) offer new 
windows and opportunities to test the proposed QCD dynamics with its pros and cons. Natu- 
rally the normal expectations run high that the bulk properties of the collision system such as 
all the momentum spectra and correlations of all produced hadrons should follow the strictures 
of QCD. But this not definite and concretely-shaped knowledge about how this actually happens 
and to what extent the process could be understood in the perturbative and non-perturative 
domains. The issues involved here still remain, to a considerable extent, quite open [2], [3]. Thus, 
having been somewhat repulsed by the so-called standard approach, we try here to explain some 
crucial aspects of measured data on pp reactions at the LHC range of energies with the help of 
some alternative approach. Our main thrust would be on the properties of time-tested familiar 
observables like transverse momenta spectra, rapidity distributions, the ratio-behaviours and av- 
erage transverse momenta (< px >) for the charged secondaries in high energy pp interactions. 
Comparison with some other model would be made whenever possible. 

The organisation of the paper is as follows: In Section 2, we provide a brief outline of the model 
chosen for study. In Section 3, the results obtained by the model-based study are presented. In 
Section 4, we end up with a discussion on the results and the observations made in Section 3 
and the conclusions. 

2 The Approach: An Outline 

This section gives a brief overview of the model-based features for the production mechanism 
of the secondary hadrons in nucleon-nucleon (p + p) interaction in the context of the Sequential 
Chain Model (SCM). According to this Sequential Chain Model (SCM), high energy hadronic 
interactions boil down, essentially, to the pion-pion interactions; as the protons are conceived in 
this model as p = (7r + 7r i?), where $ is a spectator particle needed for the dynamical generation of 
quantum numbers of the nucleons [l]-[9]. The production of pions in the present scheme occurs as 
follows: the incident energetic 7r-mesons in the structure of the projectile proton(nucleon) emits 
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a rho(g)-meson in the interacting field of the pion lying in the structure of the target proton, 
the £>-meson then emits a 7r-meson and is changed into an omega(u;)-meson, the w-meson then 
again emits a 7r-meson and is transformed once again into a g>-meson and thus the process of 
production of pion-secondaries continue in the sequential chain of q-uj-tt mesons. The two ends 
of the diagram contain the baryons exclusively [I]-[9|. 

For K + (K~ )or K°K° production the model proposes the following mechanism. One of the 
interacting -/r-mesons emits a ^-mesons; the £5-mesons in its turn emits a </>°-meson and a 7r-meson. 
The 7r-meson so produced then again emits g and (fP mesons and the process continues. The 
mesons so produced now decays into either K + K~ or K°K° pairs. The q-tt chain proceeds in 
any Fenymann diagram in a line with alternate positions, pushing the (fP mesons (as producers 
of K + K~ or K°K° pairs) on the sides. This may appear paradoxical as the production cross- 
section is generally smaller than the KK production cross-section; still the situation arises due 
to the fact that the resonances produced in the collision processes will quickly decay into KK 
pairs, for which the number of i^ will be lower than that of the KK pairs. Besides, as long as 0° 
mesons remain in the virtual state, theoretically there is no problem, for <jpK + K~ ( or (j)°K K°) 
is an observed and allowed decay mode, wherein the strangeness conservation is maintained 
with the strange-antistrange coupled production. Moreover, (fpK + K~ ( or cfpK^K®) coupling 
constant is well known and is measured by experiments with a modest degree of reliability. And 
we have made use of this measured coupling strength for our calculational purposes, whenever 
necessary. It is assumed that the K + K~ and K°K° pairs are produced in equal proportions 
|4J-[9J. The entire production process of kaon-antikaons is controlled jointly by the coupling 
constants, involving Q-ir-cj) and (fr°-K + K~ or (j) -K°K . 

Now we describe here the baryon-antibaryon production. According to the SCM mechanism, 
the decay of the pion secondaries produces baryon-antibaryon pairs in a sequential chain as 
before. The pions producing baryons-antibaryons pairs are obviously turned into the virtual 
states. And the proton-antiproton pairs are just a part of these secondary baryon-antibaryon 
pairs. In the case of baryon-antibaryon pairs it is postulated that protons-antiprotons and 
neutrons-antineutrons constitute the major bulk, Production of the strange baryons-antibaryons 
are far less due to the much smaller values of the coupling constants and due to their being much 
heavier. 

The field theoretical calculations for the average multiplicities of the tt, K and p-secondaries 
and for the inclusive cross-sections of those secondary particles deliver some expressions which 
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we would pick up from [l]-[9]- 

The inclusive cross-section of the 7r~-meson produced in the p + p collisions given by 



\ w ^-x = r 7r- exp(-2.38 < n n - > pp x) exp(- — -) , (1) 

dp 6 yy (N% ) < n n - >pp(l- x) 

Pt 

with 



< n w + > pp ^ < n w - > pp ^ < n n o > pp ^ l.ls 1/5 , (2) 

where T w - is the normalisation factor which will increase as the inelastic cross-section increases 
and it is different for different energy region and for various collisions, for example, = 90 

for Intersecting Storage Ring(ISR) energy region. The terms pr, x in equation (1) represent 
the transverse momentum, Feynman Scaling variable respectively. Moreover, by definition, 
x = 2pL/y/~s where pl is the longitudinal momentum of the particle. The s in equation (2) is 
the square of the cm. energy. 

l/p T R of the expression (1) is the 'constituent rearrangement term' arising out of the partons 
inside the proton which essentially provides a damping term in terms of a power-law in pt with 
an exponent of varying values depending on both the collision process and the specific pr-range. 
The choice of Nr would depend on the following factors: (i) the specificities of the interacting 
projectile and target, (ii) the particularities of the secondaries emitted from a specific hadronic 
or nuclear interaction and (iii) the magnitudes of the momentum transfers and of a phase 
factor (with a maximum value of unity) in the rearrangement process in any collision. And 
this is a factor for which we shall have to parameterize alongwith some physics-based points 
indicated earlier. The parametrization is to be done for two physical points, viz., the amount of 
momentum transfer and the contributions from a phase factor arising out of the rearrangement 
of the constituent partons. Collecting and combining all these, we proposed the relation to be 
given by [10] 

N R ; = 4 < N part > x /3 e , (3) 

where < N par t > denotes the average number of participating nucleons and 9 values are to 

be obtained phenomenologically from the fits to the data-points. In this context, the only 

additional physical information obtained from the observations made here is: with increase in 

the peripherality of the collisions the values of 9 gradually grow less and less, and vise versa. 

Similarly, for kaons of any specific variety ( K + , K~, K° or K° ) we have 

d 3 a. . 1 -1.33pl . 
E-j-j\ pp _, K - x = T K - exp(-6.55 < n K - > pp x) exp( ^) , (4) 



Pt 



with |r^-| = 11.22 for ISR energies and with 

< n K + > pp =< n K - > pp =< n K o > pp =< n^o > PP = 5 x l(rV /4 . (5) 

And for the antiproton production in pp scattering at high energies, the derived expression for 
inclusive cross-section is 

E-rrlpp^px = rpexp(-25.4 < n p > pp x )-—e W ( ^fjf *±) , (6) 

^ Prp <C Tip >p'p (1 - X) 

with \Tp\ = 1.87 x 10 3 and mp is the mass of the antiprotons. For ultrahigh energies 

< n p > pp ^< n p > pp ^ 2 x 10~ 2 s 1/4 . (7) 

3 The Results 

Now let us proceed to apply the chosen model to interpret some recent experimental results of 
charged hadrons production for p + p collisions at different energies. Here, the main observables 
are the inclusive cross-sections or invariant yields, rapidity distributions, ratio behaviour and 
the average transverse momenta. 

3.1 Inclusive Cross-sections 

The general form of our SCM-based transverse-momentum distributions for p+p — > C~ +X-type 
reactions can be written in the following notation: 

(i 3 cr 1 . 2 . 

E-j-j \p P ^c-x = ate- exp(-Pc- x Prh ( 8 ) 

P P T R 
The value of a n - , for example, can be calculated from the following relation: 

a n - = T n - exp(-2.38 < n n - > pp x) (9) 

The values of (a 7T -) pp , (iVJj ) pp and (/3 7T ~) PP for different energies are given in Table 1. The 
experimental data for the inclusive cross-sections versus pxiGeV/c] for ir~ production in p + p 
interactions at y / i]vlv = 62.4 GeV and 200 GeV are taken from Ref. [11] and they are plotted 
in Figs. 1(a) and 1(b) respectively. The production of 7r~, K~ and p at mid-rapidity in proton- 
proton collisions at ^snn = 900 GeV has been plotted by lines in Fig. 1(c). Data are taken from 
[12]. For the data for charged particle distribution Ed 3 N c h/dp 3 = 1 / (2ttpt) E / p(d 2 N c h / drjdpT) 
at energies ^/snn = 546 GeV and = 900 GeV we use references [3], p3]. And for LHC 



data for charged particle distribution for energies ^snn= 0.9 TeV, 2.36 Tev and 7 Tev we use 
references [14], [15]. They are plotted in Figure 2 and Figure 3 respectively. The solid lines in 
those figures depict the SCM-based plots. As the main variety of the charged particles coming 
out are the pions, we use here eqn.(l) for calculational purposes. The 'NSD'-term, used by the 
experimentalists, has the meaning of non-single diffractive collisions [16j . 

A comparison between the SCM-based results and the Tsallis parametrization is done for 
energies s/snn= 0.9 TeV, 2.36 Tev and 7 Tev. The Tsallis parametric equation [IT], [18] is 
given hereunder 

E d?N ch _ 1 E d 2 N ch _ c dN ch | ET 
dp 3 2irpT p drjdpT dy nT ' 



with y = 0.5 ln^E 1 + p z )/(E — p z )], Et = y m 2 + p\ — m and m is the charged pion mass. The 
dotted lines in Fig. 3 depict the Tsallis parametrization. 

Moreover, in Fig. 4(a) and 4(b), we have plotted theoretical values Nr and /3 versus ^Jsnn 
respectively. 

Similarly, by using eqn.(4), eqn.(5), eqn.(6) and eqn.(7), the values of (ax-)pp, {N^ ) pp , 
(Pk-)pp an d (ctp)pp, (Nr)pp> (Pp)pp are given in Table 2. The experimental data for the inclusive 
cross-sections versus pT[GeV/c] for K~ and p production in p + p interactions at ^/snn = 62.4 
GeV ,200 GeV are taken from Ref. [IE] and for y/sjfif = 900 GeV we have used Ref. [T2]. They 
are plotted in Figure 1(a), Figure 1(b) and Figure 1(c) respectively. The solid lines in those 
figures depict the SCM-based plots. 

3.2 The Rapidity Distribution 

For the calculation of the rapidity distribution we can make use of a standard relation as given 
below: 



dy J dp 



(E- n ^)dp T (11) 



In Table 3 we had made a comparison between experimentally found dn/dy for tt , K and p in 
p + p collisions for RHIC and LHC energies y^JViV = 62.4 GeV, 200 GeV and 900 GeV and the 
SCM-based calculated results. Data are taken from refs. [TT], [12] . The theoretically calculated 
results are coming out with the help of eqn. (1), eqn.(4), eqn.(6) and eqn. (11). 

Similarly, for LHC-energies, by using eqn.(l) and eqn. (11), the SCM-based dN c h/dr] will be 
given hereunder 

= 3.64exp(-0.007sinhr ? ) for y/I^ = 0.9 TeV, (12) 
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= 4.75exp(-0.009sinh7 ? ) for y^viv = 2.38 TeV, (13) 



6.28 exp(-0. Oil sinhr?) for = 7 TeV. (14) 



and 

dNch 
drj 

In Fig. 5 we have plotted dN c h/drj vs. rj at three LHC-energies ^/snn = 0.9 TeV, 2.38 TeV 
and 7 TeV. The reconstructed data points for Fig. 5 are from Refs. |14| . |16j . |19j . Lines in the 
Figure are the outcomes of eqn. (12), eqn.(13) and eqn.(14) respectively. 

3.3 The Ratio-behaviours for Different Secondaries 

The nature of the relation of K/tt ratios with the SCM, presented in the previous work [2D] , 
would be written in the following form 

— = 5.4x 10- 2 (v^) ai . (15) 

TT 

Fig. 6(a) shows the nature of rise of K/tt ratio in the light of SCM-based above relation (eqn. 
(15)). Data are taken from Ref. [21] . 

Similarly, in Fig. 6(b), we have presented the p/ir ratio for the RHIC and LHC-data.[21J The 
SCM-based calculations are done on the basis of eqn.(l) and Table 1. 

3.4 < pr > Values 

Next we attempt at deriving model-based expression for < px >■ 

The definition for average transverse momentum < px > is given below. 

< p T > C = PT(m ; n) W r , (16) 

r PT (max) „ d 3 a ° j 2 
Jpr(min) dp 3 "T 

The line in the Fig. 7 depicts the SCM-based calculated results of the average transverse mo- 
mentum < px > versus the cm. energy -^/snn- The theoretical calculations are done on the 
basis of uses of eqn.(l), Table 1. Data are taken from [3], [T4], [22], [23]. 

4 Discussions and Conclusions 

Let us make some general observations and specific comments on the results arrived at and 
shown by the diagrams on a case-to-case basis. 
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a) The measures of inclusive cross-sections against transverse momenta (jpr) obtained on the 
basis of the SCM for the pions, kaons and protons for the RHIC energies ^snn = 62.4, 200 
GeV and for LHC energy ^Jsnn = 900 GeV are depicted in Fig. 1. They describe a modest 
degree of success. 

Of the various types particles produced, the 7r-mesons, constitute, the near totality of the 
secondaries. So, in calculating the charged hadrons yields for different transverse momenta, on 
the basis of the model, for different energies ranging from ^snn = 546 GeV to 7 TeV, we use 
eqn. (1). The results are shown in the Fig. 2 and Fig. 3. Moreover, we have compared the 
Model-based results with Tsallis parametrization in Fig. 3. The outcomes of these plots are fairly 
satisfactory. 

b) The theoretically calculated Nr and /3-factors of eqn. (8) for different ^/sjvjv's have been 
plotted in Fig. 4. With the inclusion of the power law form arising due to the physics of partonic 
rearrangement factor, the model has turned effectively into a mixed model. And these two 
factors corroborate effectively the "soft" and "hard" regimes without any extra effort. 

c) The calculations of rapidity distributions for pions, kaons and protons at different RHIC 
and LHC energies, on the basis of the model, have been done with the help of eqn. (11). The 
calculated values are compared with the experimental ones and they are shown in Table 3. The 
theoretical values are in fair agreement with the experimental data. 

Similarly, in Fig. 5, we have plotted the pseudorapidity distributions for different LHC ener- 
gies. Here, the model modestly reproduce data consistently. 

d) The agreements between the measured data on K/tt and p/tr ratio and the theoretical 
SCM plots (shown in Fig. 6) for different energies are strikingly encouraging. 

e) Fig. 7 shows the plots of < px > vs. ^snn. The initial indication of the SCM-based 
theoretical plot shows a modest agreement with the data. 

Finally, we conclude from the analysis of the results given above with the following statements: 
The model applied here gives fair descriptions of the pr-spectra of all the light secondaries 
or charged hadron. Some disagreements are observed at pseudorapidity distributions at LHC 
energies. The model might require some finer adjustments to cope with the data for very high 
energy nuclear collisions. However, on an overall basis, our model is in fair agreement with the 
latest pp-collision results obtained from the uptodate LHC experiments. This factor is really 
of high important to us. Furthermore, the observables or the physical aspects that we have 
reckoned herewith form a clear continuum from the old ISR experiments to the recent Large 
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Hadron Collider studies via the intermediary BNL-RHIC results cropped up over the 1st decade 
of this century. In so far as the rolls of the other models (including QCD versions) are considered, 
the results do neither speak very high about them; rather they cast doubts on the suitability of 
them in applying at this LHC energy band. And this certainly spurs us on to take up further 
studies on the SCM proposed earlier and applied in the present study. 
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Table 1: Values of (a 7T -) pp , (iVjj ) pp and (/3 7r -) pp for it productions in p + p collisions at 
^/IWn=Q2A, 200, 546, 900, 2380 and 7000 GeV 



y/SNN 


{ a TT- )pp 






62.4 GeV 


0.545 


3.327 


0.468 


200 GeV 


0.907 


3.867 


0.293 


546 GeV 


0.153 


3.931 


0.172 


900 GeV 


0.135 


4.154 


0.128 


2380 GeV 


0.166 


4.235 


0.085 


7000 GeV 


0.355 


4.366 


0.075 



Table 2: Values of (a K -) pp , (N% ) pp , (f3 K -) pp and (a p ) pp , {N P R ) PP and (f3 p ) pp for K , p pro- 
ductions in p + p collisions at ^5^^=62.4, 200 GeV and 900 GeV 





(a K - ) pp 




{Pk- )pp 


62.4 GeV 


0.245 


2.527 


0.591 


200 GeV 


0.235 


3.017 


0.417 


900 GeV 


0.047 


1.544 


0.248 


\j $NN 


i a P ) PP 


(K) PP 


{Pp) PP 


62.4 GeV 


0.215 


1.527 


0.618 


200 GeV 


0.115 


2.117 


0.426 


900 GeV 


0.027 


1.344 


0.248 



Table 3: Comparisons of experimental dn/dy with the SCM-based theoretical ones for it , K , 
p productions in p + p collisions at ^/satat=62.4 GeV, 200 GeV and 900 GeV 





particle 


dn/dy 
(Experimental) 


dn/dy 
(Theoretical) 


62.4 GeV 




0.900±0.063 


0.840 




K~ 


0.103±0.005 


0.100 




p 


0.037±0.003 


0.035 


200 GeV 




0.824±0.053 


0.841 




K~ 


0.067±0.003 


0.065 




p 


0.022±0.002 


0.018 


900 GeV 




1.485 ± 0.004 


1.241 




K~ 


0.182 ± 0.004 


0.175 




p 


0.079 ± 0.002 


0.068 
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0.001 r 



0.0001 




1.5 2 2.5 

p T (GeV/c) 



(c) 



Figure 1: Plots for ir, K and proton-production in p + p collisions at RHIC energies (a) ^/snn = 62.4 
GeV, (b)^/sWN = 200 GeV and (c) ^/sWn = 900 GeV. Data are taken from Ref. [TT] for Figs, (a) and 
(b), while for Fig. (c) from Ref. [12] . Solid lines in the Figures show the SCM-based plots. 
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10000 




1e-006 1 1 ' 1 ' 1 1 

1 2 3 4 5 6 

p T (GeV/c) 

Figure 2: Transverse momentum distributions 
of the invariant cross section of charged parti- 
cles in NSD p + p collisions at \y\ < 2.5 for 
energies ^snn = 546 GeV and 900 GeV. Ex- 
perimental data are taken from [3], |13j . Lines 
show the theoretical plots. 



p+p = h + "+X (LHC-CMS Data) 




1e-007 1 1 1 1 1 1 1 

1 2 3 4 5 6 

p T (GeV/c) 



Figure 3: Charged-hadron yield for energies 
= 0.9 TeV , 2.38 TeV and 7 TeV in 
the range r\ < 2.4 in NSD events as a function 
of pt', Data are taken from CMS collaboration 
Q3]-[Ti]. Solid lines in the Fig. represent SCM- 
based results while the dashed lines show Tsal- 
lis fit. 
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Figure 5: Plot of dN ch /dr] vs. rj for p + p collisions at ^/sWn = 0.9 TeV, 2.38 TeV and 7 TeV. The 
reconstructed data points are from Refs. [14], [16], [19]. The lines in the figure depict the theoretical 
results for different energies. 
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Figure 6: Presentation of plots of (a) K/ir and (b) p/n at different center-of-mass energies. Line and 
rhombus represent the SCM-based results in (a) and (b) respectively against the data sets taken from 
Ref. [UJ 
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Figure 7: Plot of average transverse momentum as function of y/s. Data are taken from [3], [T4], [2"2"] . 
[23] , Line shows the SCM-based calculations. 
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